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In order to calculate cross sections with a large number of particles/jets in the final state at next-to- leading 
order, one has to reduce the occurring scalar and tensor one-loop integrals to a small set of known integrals. In 
massless theories, this reduction procedure is complicated by the presence of infrared divergences. Working in 
n — 4 — 2e dimensions, it will be outlined how to achieve such a reduction for diagrams with an arbitrary number 
of external legs. As a result, any integral with more than four propagators and generic 4-dimensional external 
momenta can be reduced to box integrals. 



1. Motivation 

The study of processes with multi-particle/jet 
final states is of major importance not only to test 
perturbative QCD, but also to obtain a precise es- 
timate of the background for search experiments. 
At the Tevatron, multi-jet cross sections with up 
to six jets have been measured and compared to 
tree level QCD predictions [ ||. However, leading 
order predictions are very unstable with respect 
to variations of the renormalization and factor- 
ization scales, such that the NLO corrections are 
necessary to obtain reliable predictions. 

The calculation of e.g. a 2 — > N — 2 parton pro- 
cess at NLO requires the knowledge of iV-point 
scalar and tensor one-loop integrals with massless 
propagators and external legs, for e + e~ — > N jets 
one needs (N + l)-point integrals with one off- 
shell external leg. These integrals have to be fur- 
ther processed by reducing the tensor integrals 
to scalar integrals and, for N > 5, reducing the 
scalar iV-point integrals to known integrals with 
less propagators. 

Of course, reduction methods have been 
worked out and successfully applied before [ 0|— [ 
but mostly for integrals with massive prop- 
agators, as needed e.g. for electroweak radiative 



corrections. In contrast, when dealing with mass- 
less partons, the presence of infrared divergences 
requires a different approach. For this case, re- 



duction methods have been worked out [ g-[ 14 
which are valid for up to 6 external legs. How 
to extend these results to more than 6 external 
legs was not clear until recently [ || due to prob- 
lems stemming from the inversion of kincmatical 
matrices which were singular if the external mo- 
menta were kept in four dimensions. On the other 
hand, the use of helicity techniques requires 4- 
dimensional external momenta. The main ideas 
how to overcome this problem will be outlined in 
the following, leading to reduction formulas for 
scalar and tensor integrals valid for an arbitrary 
number of external legs. 

2. Reduction of scalar integrals 

Consider the scalar integral 
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in n = 4 — 2e dimensions with TV massless propa- 
gators, depicted in Fig. [j]. 

The corresponding expression in Feynman pa- 
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Figure 1. The one-loop N-point graph. 



rameter space reads 
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The kinematic information is contained in the 
matrix S which is related to the Gram matrix 
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The external momenta are kept in four dimen- 
sions; momentum conservation implies = 0. 
The integral 1^ can be split into a singular part 
with less propagators and a finite, but formally 
higher dimensional integral as 
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where the coefficients bi have to fulfill the follow- 
ing equations: 
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The integral Idiv in (0) is a sum over reduced in- 
tegrals where one propagator could be cancelled. 
Our aim is to show that the higher dimensional 



integrals I N +2 in (|3|) , which are very hard to cal- 
culate for N > 4 and hinder a reduction, drop out 
for TV > 4. To this end we have to solve Eqs. (Q) 
for the parameters bi. Note that for N = 5, ex- 
pression ([|) is of order e, which means that Ifi n 
can be dropped right away in that case because 
(n + 2)-dimensional integrals are finite for N > 4. 
Since only four 4-dimensional external vectors 
can be linearly independent in Minkowski 
space, rank(S')=min(A r ,6) and rank(G)=min(./V- 
1,4). Hence the case N = 6 is special since S 
is still invertible for N = 6, whereas G is not. 
Therefore Eqs.(Q) can still be solved by inverting 
S for N = 6, whereas for N > 6 both, G and 
S have a vanishing determinant. In this case, 
the solution to (Q) can be constructed in terms 
of the so-called pseudoinverse H to G. It is de- 
fined by the properties HGH = H, GHG = G. 
Given a pseudoinverse H a solution to Gx = y 
exists if and only if y — GHy. Introducing a ba- 
sis Ei_ x 4 of Minkowski space and defining the 



matrices R and G by 
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the uniquely definedn pseudoinverse is given by 



H = R A {RR )-' G (RR 1 )" i R . 
Now Eq. ([|) can be solved for general N [ 



(G) 



2 In general H is not unique. In our case, the uniqueness 
follows from the symmetry of G. 
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Doing so one finds 

N 



J2 b i = ° for AT > 6 



such that the contribution (H) completely van- 
ishes for all N > 6. This means that all scalar 
N— point integrals, for arbitrary N, can - by re- 
cursion - be expressed in terms of box integrals. 
Explicitly, one finds the following reduction for- 
mulas: 

For N < 6, i.e. det(S) ^ 0: 



1 N 

o E ^hl 1 Iff -1, 
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2 Jv det(5) w 
For iV > 6, i.e. det(S) = 0: 

AT-1 JV-6 



(7) 



(fe) 



/=1 



jV ± N-l,ll 



k=l 



(8) 



The reduced integrals /jv— 1,1 are defined as 



J JV-1,. 
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and if = 1jv-i -HG, u; = rf . 

The (iV — 6) /3ft in (||) are free parameters 
to construct linear combinations of the vectors 
{ [/(*)} (fc = 1, . . . , N - 6) which, together with 
jj( N ~ 5 \ form a basis for the kernel of G. We 
chose 

UiN-5) =k-v/{v-K-v) 
parallel to v and the others orthogonal, 

v , u{ k=l,...,N-6) =Q 

Obviously one can choose the /3's such that (N — 
6) 6's from the set {bi, . . . , &at_i} are zero. Doing 
so one observes that can be expressed by only 
6 (N — l)-point graphs for arbitrary N > 6. This 
is the generalization to massless kinematics of a 
similar result in the infrared finite case [ |(| which 
has been derived for integer dimensions using 4- 
dimensional Schoutcn identities. 



3. Reduction of tensor integrals 

Tensor integrals can also be split into a part 
containing reduced integrals (K^) and a part 
containing higher dimensional objects (J N ')- Us- 
ing momentum conservation (rjv = 0) one can 
write: 
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are differences of reduced integrals and 



I N +2l (jl,.--,3L-2l) 



are higher dimensional integrals with Feynman 
parameters zj in the numerator. The Lorentz ten- 
sors T, Q present in the above expressions can be 
written in terms of the objects (vi = rf) 



r» -H ■ 



W" = ■ v 
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Note that the pseudoinverse H to the Gram ma- 
trix G also appears here. The important feature 
is that one always has 

G ~ 9^/2 -i* -H-r v . 

If now the Lorentz vectors (j = 1, . . . , N — 1) 
span Minkowski space, which is generically the 
case for N > 5, one immediately derives, using 
the explicit expresssion of H given above, that 



4-di: 
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This means 

-r" -H-r v = 0(e) 

and thus the higher dimensional objects in 
Eq. (^) can be dropped for N > 5. Consequently 
it follows from Eq. (||) that for the reduction of an 
TV-point integral only a small number of higher 
dimensional integrals is needed, i.e. one has to 
know J^- tlL only for 2 < L < N = 2, 3, 4. These 
can be calculated once and forever, the explicit 
expressions are 

^=| 3 ,4 - - (5 M1M2 /2 - n^)i n N +2 
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The indexed bracket is an abbreviation for the 
sum of all distinguishable distributions of Lorentz 
indices to the objects inside the bracket, e.g. 

•O1A12M3} 



XYY 



-X^Y^Y^ 1 + X^ 3 Y^ 1 Y^ 2 



For the n-dimensional part, K^ 1 '" fJ ' L , a recursion 
relation holds. Skipping all the details (see [|^]) 
one finds 
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Putting everything together, the final recursion 
formula for A-point tensor integrals reads 
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This form is not immediately suited for iteration, 
as some reduced integrals do not contain the triv- 
ial propagator k 2 anymore. Therefore one has to 
shift the loop momenta of the respective integrals 
by k — y k + 77 before iteration, leading to 
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The argument vectors are 

i? = (ri,...,rjv) 
%] = (ri, . . . ,f fc , . . . ,rAr) 
R [l] = {n+i-ri,ri + 2-n,. 



(13) 



,riv-2+i - n,0) 



where f means that the respective vector is miss- 
ing. The vector indices are understood to be 
taken cyclically symmetric with periodicity N, 
i.e. v = ri for I G {1, . . . , N}. As we as- 
sume r/v = 0, the integrals on the right-hand side 



of (13) have again at least one trivial propagator 
and thus are suited for a further reduction step. 
As a simple example, consider the rank two 3- 
point integral I(^ 1>12 (R). After the first reduction 
step and subsequent shift, one obtains 

jMirt = 1 JM2 + W /i 2 jMi 
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The recursive structure together with the shift 
can be implemented in an algebraic manipulation 
program in a straightforward manner. The recur- 
sion stops if all tensor integrals are transformed 
into expressions containing scalar integrals and 
Lorentz tensors composed of the objects given in 
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Eq. (|Tc|). Starting with an TV-point tensor inte- 
gral one finds a linear combination of scalar inte- 
grals with K < N legs. To these scalar integrals 
the reduction formulas (|?]),(|j) can be applied. In 
this way our reduction formalism solves the prob- 
lem of calculating massless iV-point Feynman dia- 
grams for arbitrary N. The generalization to the 
case where some of the propagators are massive 
will work similarly. 

4. Summary/Outlook 

In conclusion, massless scalar and tensor one- 
loop integrals for graphs with an arbitrary num- 
ber of external legs (N > 4) can recursively be 
reduced to linear combinations of (known) box 
integrals. To deal with the problem of singu- 
lar Gram matrices we introduced the concept of 
the pseudoinverse, which renders the whole for- 
malism mathematically well defined for arbitrary 
N. The recursive structure can easily be im- 
plemented in algebraic manipulation programs. 
This allows in principle for the NLO calculation 
of multi-jet cross sections with a large number of 
final states. Applications relevant for the near fu- 
ture are e.g. the NLO calculations for pp — » bb bb 
or pp — > 4 jets/photons. The most complicated 
analytical structure in such a calculation is the 
scalar 6-point function since, as shown above, no 
dangerous higher-dimensional objects have to be 
computed. With the given methods, an explicit 
analytical expression for the scalar 6-point func- 
tion has been derived and presented in [ |2| . 
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